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ABSTRACT 


„ y in initial crack size and fatigue crack growth is incorporated in a probabilistic 

model that is used to predict the fatigue lives for unstiffened aluminum alloy panels containing 
multiple site damage (MSD). The uncertainty of the damage in the MSD ZSKShv 

sizes tn The < ^Iwh^Tnfa^ “f ana ¥ ic fy drived from equivalent initial flaw 

Zlll Yfnability.in fatigue crack growth rate is characterized by stochastic descriptions of 
crack growth parameters for a modified Paris crack growth law. A Monte-CarL simulation 

and then^ows M S D h Y ™ ndorid y. selecting values from the stochastic variables 
i M ^ D cr ? cks with a deterministic fatigue model until the panel fails. 

Different simulations investigate the influences of the fatigue variability on the distributions of 
remammg fatigue lives. Six cases that consider fixed and variable conditions of initial crack size 
and fatigue crack growth rate are examined. The crack size distribution Lhibked a dorSnant 
°? f th f rem / un !f lg fatigue life distribution, and the variable crack growth rate exhibited a 

J® 45 th f , dlstnbu . tlon - to addition, the probabilistic model predicted that only a small 

percentage of the life remains after a lead crack develops in the MSD panel. y 


INTRODUCTION 


Multiple site damage (MSD) refers to the accumulation of widespread fatigue cracking in 

M s pca nb ee S pecialIydamagingwhencLks' ' 


* . r . „ % — - V uuuiu 6 m 5 iiiiuuie ai several 

adjacent fastener holes in a mechanically fastened joint. Over time, the effects of the fatigue 

amage will reduce the residual strength of an aircraft joint and may result in a loss of damage 
tolerance capability. To reliably predict the capability of an MSD joint, the influences of § 
widespread fatigue cracking and the associated crack growth mechanisms are needed. 

The development of a method that incorporates the randomness in widespread fatigue 
cracking lends itself to a probabilistic approach. The initial crack sizes in the MSD joint can be 
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represented by randomly selecting crack lengths from an analytically or experimentally derived 
distribution. The distribution of fatigue crack sizes can be estimated using an equivalent initial 
flaw size (EIFS) method proposed by Manning and Yang [1-5]. In addition to different MSD 
crack lengths, each crack is expected to grow at a different rate as observed by Virkler et al. [6]. 
Ostergaard and Hillberry [7] derived the distributions for the fatigue crack growth parameters 
that described the variability observed by Virkler. From the development of accurate 
descriptions of the crack sizes and crack growth rate as described above, a probabilistic model 
can be implemented to determine the influence of these random variables on fatigue life 
predictions. This approach can then be applied to an unstiffened panel containing MSD cracks to 
provide a better understanding of the effects of widespread fatigue cracking. A model of an 
unstiffened MSD panel was developed by Moukawsher et al. [8-1 1]. Moukawsher's model 
deterministically calculates fatigue crack growth, incorporates crack interaction effects, and 
predicts panel failure. This MSD fatigue model is integrated into a probabilistic model to 
determine a distribution of fatigue lives for a panel containing MSD cracks. 

The objectives of this study are: (1) to characterize the variables associated with widespread 
fatigue cracking, (2) to develop a probabilistic method that incorporates the randomness for these 
variables, and (3) to evaluate the influence of fatigue variability on life predictions for an 
aluminum alloy panel containing MSD cracks. The present goal is to observe trends in fatigue 
lives of panels that contain both multiple site damage on the scale observed in aging aircraft and 
realistic crack growth properties of the material that was used to build these aircraft. 


MATERIALS AND METHODS 


Multiple Site Damage Model 


The multiple site damage model considered here is a panel 787 mm wide by 2.3 mm thick 
containing 30 holes 4.76 mm in diameter (Figure 1). The holes are equally spaced and aligned 
perpendicular to the applied loading. This geometry was also considered by Buhler et al. [11] in 
a parametric study which employed the deterministic MSD analysis developed by Moukawsher 
et al. [8-10]. For this study the panel material considered was 2024-T3 aluminum alloy. The 
initial damage state was through-the-thickness cracks emanating from each side of the 30 holes. 
The panel was subjected to a constant amplitude remote tension of 68.9 MPa and a stress ratio of 
0.1. Failure of the panel was considered to be fracture. A residual strength criterion based on 
Swift's ligament yield model [12] was used to predict failure. Moukawsher et al. [9,10] 
experimentally tested 2024-T3 aluminum alloy panels with MSD and observed good agreement 
with the residual strength criterion. 


Variability in Initial Crack Size 


The initial fatigue cracks at the holes in the panel were represented by an analytically derived 
distribution of crack sizes. The Manning and Yang approach for predicting fatigue cracking in 
aircraft from equivalent initial flaw size distributions was followed [1-5]. An EIFS distribution 
represents a population of artificial cracks, which when grown forward will be similar to 
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^=Qa(t) 


( 1 ) 


^tlf!^ d fr t)/dt . iS th l Cfack growth rate ’ a is *e crack size in millimeters, Q is the crack growth 

^ ' ime ' ^ EIFS ™lues are & 


F( a ) = exp J - [ I n ( x u^ a ) j a 


( 2 ) 


where x u is an upper-bound, a is the scale parameter, and <$> is the shane parameter. 

Widespread fatigue cracking in aircraft is predicted by growing the EIFS distribution 

gS ?u Se f ?^f f rn servi f time t=T T° « ^ 

t SCGC ) for a given stress region of an aircraft is estimated with the analytical 
i S row th rate equation (equation 1). A reasonable approximation is achieved bv fitting the 

s^ y nSlovl- atl0n m tW0 ' SegmentS l ° 016 SCGC ac « to a refe?en^cmTsiS £)£ 


^?=Qi ’ abl 

a(t) < ao 

(3) 

“jf = Q2 ' a b2 

a(t) >= ao 

(4) 


Tire crack growth coefficients Q, and Qj for a given SCGC are estimated from equations 5 and 


Ql=lta 


^1 


Q2 = §2 ' cr 


V 2 


(5) 


( 6 ) 


where a is the maximum stress level of the SCGC in MPa, and, 5., s 2 % are emnirical 

3 S u an S based °, n crack Slze data observed fractographically. The crack’ growth exponents bi 
and b 2 were equal to one m this study as suggested by Manning and Ynn? n OnL th*> w 

anaWt^ C p efflCientS f haVe e ” estimated for a given SCGC, the EIFS distribution can be ° 
analytically grown forward to determine the distribution of crack sizes at some service time t=T. 

low^ wf^ldn P m aC ^ a p!l?? d H^ Y u n ? and Manning in the dura bility analysis of a fighter 
lower wing skin [1] The EIFS distribution parameters and empirical constants were determine 

tSS ^'°,f? PhlC data 5615 f ? r lab0rat0I y s P“imens under spe^mSSg S reluS 
EIFS distnbution parameters for equation 2 were x u =0.762 mm, 0=1.716, <4 308 and the 

constants for equations 5 and 6 were |i=2.227E-19, | 2 =6.288E-8, 4> 1= 6.374, and 4> 2 =1 546 The 

wW i r -°T h coefflcients (Ql and Q 2 ) were then estimated for ten different stress regions in the 
wmg skm from equations 5 and 6. The analytically derived EIFS distribution can thin be grown 
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forward by integrating equations 3 and 4 for each wing skin region to determine a distribution of 
crack sizes at any service time t=T. 

Data reported by Yang and Manning [1] for the fighter lower wing skin were used in this 
study to generate the prediction for the extent of widespread fatigue cracking. The EIFS 
distribution was grown forward to 10,000 flight hours for stress region VI. The resulting fatigue 
crack sizes were then fit to a Weibull compatible distribution function (equation 2) with 
parameters x u =4.27 mm, a=1.995, and <(>=6.921. This analytically derived distribution is used to 
represent the initial crack sizes in the MSD panel. 


Variability in Fatigue Crack Propagation 


The variability in fatigue crack growth behavior was observed by Virkler et al. [6] for 68 
replicate, constant amplitude tests of 2024-T3 aluminum alloy panels (Figure 2). Virkler 
measured the fatigue crack growth (a vs. N ) in center cracked panels cut from the same 
material lot. The data were readily described by a modified Paris law [7] 

^ = C-(AK-AK th ) m O) 

where the crack growth rate parameters C and m are material constants, da/dN is the crack 
growth rate in inches/cycle, AK is the stress intensity range in ksiVni, and AK t h is the threshold 
stress intensity. The threshold stress intensity value used by Ostergaard and Hillberry [7] and in 
this study was 2.75 ksiVni. Ostergaard and Hillberry [7] used a finite integral optimization (FIO) 
method to determine the fatigue crack growth rate parameters such that upon integration of 
equation 7 the original a vs. N data were accurately reproduced. The FIO method optimizes 
each parameter set by determining the Cj and mi that give the best-fit to the a vs. N data 
according to a prescribed minimum error function. 

From the 68 sets of optimized values for the crack growth parameters Cj and mj, Ostergaard 
used a least-squares regression analysis of the log m vs. log C data to show that the mean crack 
growth behavior could be described by equation 8, 

log Q = bo + bi • log mi (8) 

with mean crack growth parameter Q. The least-squares estimates were bo = -5.7792 and bi = 

-4.6150. The deviations of the fitted value log Q (from equation 8) compared to log C were used 
to describe the variability in crack growth rate. This variability was characterized by the 
parameter F defined below 



(9) 


log Fj = - (log Q - log Cj ) 


( 10 ) 
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resulfs arc r f pres ? n,s «ie least-squares residuals from the regression analysis. The 

fatigue cmck growtim-afe nt “** gr ° W ' h Parame ' erS mi and F ‘ tha ‘ describe ,he variability in the 

3-p^n^^ " da ' a • b est-fi, of the F parameter to the 


f<F) = 


: ^feexp(-i f .„ g (F-r)-a )2) 


( 11 ) 


s hape and threshold parameters were estimated to be a = -0 717 B = 0 075 and 
«WoS2f' ” parameter ShOWed a best - fit '° ,he 2 -P^meter loAo^ai 5 ' 


f(m) = 


exp (~ ITT ^°g m - oc } 2 


\ 2 P 


) 


( 12 ) 


where the scale and shape parameters are a = 0.35508, and p = 7. 1 1 62E-4 respectively The 
£uuu£ ta S5SSi*SST^ 2 " USed “ the probabilistic » -preset n, the 


Probabilistic Modeling 


the^itiSSe SocL a s SP ^^?nf^ liStiC m ?, d ? 1 j ng \ s to use a mecha ™stic approach that describes 
g e process. The MSD fatigue model developed by Moukawsher et al T8-1 11 nrovide<; n 

MSD Sacks ‘ ffST ‘° predic, ; he , fati ^ e life of a rcmote 'y loaded pand containing 

effects clkuSrionnf th?? St t e locat i on ’ len § th ’ f yp e of crack, and crack interaction 
eirects in the calculation of the stress intensity factor (aK) solution for each crack The 

corresponding crack growth rate is calculated from the da/dN vs aK relationship in this case a 

modified Pans law. The remaining fatigue life of the MSD panel is predated by’ 

criterSn. 1StlCally §r0Wing the cracks untl1 the P anel faiIs according to a residual^strength 

j A .^onte-Carlo simulation incorporates the fatigue variability with the MSD fatigue model 
desenbed above to determine a distribution of remaining fatigue lives for the MSD panel This 
probabilistK approach is used ,o investigate the influences of the variaWHtyln fttlgL Chador 

considered 6 wer^ Tf mea c ° m P ann f dlfferent simulation conditions. Thecrack conditions 
ed ?'• f an size cracks emanating from both sides of each fastener hole ?) a 
random symmetric crack size emanating from both sides of each fastener hole and 3) a random 

3S““ k from both sides of each fastened Me Thecmck gr3 mTe 

conditions considered were: 1) mean crack growth parameters 2) random crack orowth 

P“ rS f ° r t f ach pa ^ el (homogeneous panel), and 3) random crack growth parameters at 
every ligament for each panel (heterogeneous panel). parameters at 

™ v Monte-Carlo simulations of 1000 trials each were investigated for the three crack size 
conditions and the three fatigue crack growth rate conditions. The influence of the fatigue 
variability was considered by evaluating simulations with different conditions of one fatigue 
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variability, while holding the other variability the same. The different simulation cases are 
described below, and are identified in Table 1. 

• Case 1 : mean crack sizes emanating from both sides of each hole and a randomly 
selected crack growth rate for each panel. 

• Case 2: randomly selected symmetric crack sizes and a mean crack growth rate for 
each panel. 

• Case 3: randomly selected unsymmetric crack sizes and a mean crack growth rate for 
each panel. 

• Case 4: randomly selected symmetric crack sizes and a randomly selected crack 
growth rate for each panel. 

• Case 5: randomly selected unsymmetric crack sizes and a randomly selected crack 
growth rate for each panel. 

• Case 6: randomly selected symmetric crack sizes and a randomly selected crack 
growth rate at every ligament for each panel. 

By direct comparison of the fatigue life distributions for these cases, the fatigue life variability 
can be quantified to provide a better understanding of the uncertainties in the MSD panel. 


Table 1. Monte-Carlo Combinations of the MSD Panel 


Case 

Crack size 
parameter 

Material 

parameter 

1 

mean crack size 
every hole 

variable crack growth 
rate each panel 

2 

variable symmetric 
crack size every hole 

mean crack growth 
rate each panel 

3 

variable unsymmetric 
crack size every hole 

mean crack growth 
rate each panel 

4 

variable symmetric 
crack size every hole 

variable crack growth 
rate each panel 

5 

variable unsymmetric 
crack size every hole 

variable crack growth 
rate each panel 

6 

variable symmetric 
crack size every hole 

variable crack growth 
rate every hole 
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RESULTS AND DISCUSSION 


comnarpuh ^^ deviations are plotted m a bar chart for each case (Figure 9) Case 1 vs 4 

rate condition crack sizes and tbe randomly selected crack sizes for the same crack growth 
sdLted crS; SI i comparison, the mean fatigue life decreased significantly for the randomly 
hot? Ac k S Cracks larger than the mean value interacting with the adjoinin" crackeS 
S £?. a resub ’ a dominant crack would link-up to form a lead Srack, SSKd 

fail soon thereafter as a result of the loss in residual strength. Case 2 vv 4 coirmares the 
Se ioinv gr ° M rate .™ th the randomly selected crack growth rate for the same inS craS 
a .^ XOa ‘ .^° S1 §mficant difference in the mean fatigue lives occurred - however ihe 
“ creased with the addition of the variable crack growfh ’ Thfs su-Sts that 

of the dStriSonSs iffeSed n S 1VC S the . variabilit y in fati gue crack growth rate, bSTthe shape 
- ts affected. Case 4 vs. 6 compares the variable crack growth rate for each 1 

aad the vanable crack growth rate for each ligament in the panel Both cases exhibited 

co^mrefth lbUtl ° nS WhlCh su PP ort the previous observations. Case 2 vs. 3 and Case 4 vs 5 
compares the symmetric crack sizes and the unsymmetric crack sizes These comparisons Sow 
an increase m the mean fatigue lives and a decrease in the standard deviation fonhe 
„ vm y™? etnc crack conditions. The increase in mean life is expected since a lar^e crack that is 

unsymmetric crack condition. Overall the crack 
3 d °ff inant ef l ec l on the fati § ue life predictions, and the^We^k 

"“”ber of holes and the Jigue v^SS^ life 

resid^'stre^gA^f the^anel^ < ^is^havio^asinvestigat^ a by c^culldfngThe p^cent°^ tbe 

thmurflTa 1 r after the fi u st I,nk '“p- Histogram plots of the percentages areshown in Figures 10 
YSPlSji 4 I Comparing the means and standard deviations for each case in the bar chart '(Figure 
15) shows that approximately 5% of the life remains after the development of' a lSd crack ¥hfc 
suggests that the majority of the life is spent propagating the initial cracks' therefore a reliable 
analysts must include a good representation of the input'variables. The uniform crac'kina n Cure 
1 never developed a lead crack and the panel failed from net section yielding 8 

, h ' Tjs utlIlty ‘. n analyzing trends between different cases depends on appropriately reoresentino 
ue .^i lines m the MSD panel. It was determined that for the Sved 

thcmohif- d r bUt,0n ’i hat a PP roxlma tely 33% of the cracks would propagate y Thisls a result of 
n Pans crack growth equation (equation 7) preventing cracks from growing afstress 

10 (S flf«foh W * e threSh ° W Val “ e - The derived cra = k size diftribution to”h! MSD pldlt 

traW? fllghtboars was compared to observed crack size data reported by Lincoln r 13 /for 19 

^ealinT2fio%ofm n , ?? d " Stra ? ive tear down and ias P ecd «n oSner wings 

a .JT? 25 ; 50 ^. ofth c fastener holes had detectable cracks. Additionally the initial crack size 

?Zf t butl ° n f° r . thls study compares favorably with the fitted distribution to the observed cracks 

we U S 16 > distributions of fatigue 

aTfoeSS^nthe^ a!“eef ^ “ 1976 W “ Ch iS 0f the " a “ mintage 

Although the cuirent paper has expressed the results of the probabilistic MSD analvsis in 
terms of remaining fatigue life distributions, the approach could also be used to determine the 
variation in residual strength as a function of elapsed cycles. In that case, Monte-Carlo 
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simulations would be performed as before to predict distributions in MSD crack patterns at 
various predetermined life intervals. The various crack patterns would then be combined with a 
failure criterion to compute the distribution in residual strength for a given point in lire. 
Determining the distributions in residual strength as a function of elapsed cycles would, then, 
give an indication of how the safe load carrying capability of the structure decays with service 

usage. 


SUMMARY 


A probabilistic model that includes the variability in crack size and fatigue crack growth rate 
has been developed for investigating the fatigue lives of MSD panels. These variables were 
defined from examples that represent MSD in aging aircraft and represent the properties of the 
materials that was used to build these aircraft. Six Monte-Carlo simulations were run with 
variations of fixed and variable conditions for different random variable combinations. The 
MSD crack size distribution had a dominant effect on the fatigue lives, and the variability in 
crack growth rate had a lesser effect. The model predicted that the majority of hfe was spen 
propagating the MSD cracks until adjacent crack-tips link-up, and then approximately 5% of the 
panel life remains. The probabilistic model provides an improved understanding of the 
influences of widespread fatigue cracking on remaining fatigue lives. 
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Symmetric Cracks 


Figure 1: 


Figure 2: 



Linked Cracks 


Material Properties: 

Aluminum 2024-T3 
Thickness: 2.286 mm 
Width: 787.4 mm 
Hole Diameter: 4.76 mm 
Hole Pitch: 25.4 mm 


Hill 

Applied Stress 


MSD panel. 



Virkler's data from 68 replicate fatigue crack growth tests in 2024-T3 aluminum alloy 
specimens. 
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Figure 3: Remaining fatigue life histogram. Case 1 : mean size cracks; variable growth 
parameters for each panel. 6 
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Figure 4: Remaining fatigue life histogram. Case 2: variable symmetric crack-sizes; mean crack 
growth parameters. 
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Figure 5: Remaining fatigue life histogram. Case 3: variable unsymmetric crack sizes, mean 
crack growth parameters for each panel. 
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Figure 6: Remaining fatigue life histogram, Case 4: variable symmetric crack sizes; variable 
crack growth parameters for each panel. 
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Figure 8. ^™m m g fatigue life histogram, Case 6: variable symmetric crack sizes; variable 
ack growth parameters at each hole for each panel (heterogeneous panel). 
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Figure 9: Remaining fatigue life bar chart of mean and standard deviation for Cases 1-6. 



Figure 10: Percent remaining fatigue life histogram after developing a lead crack, Case 2: 
variable symmetric crack sizes; mean crack growth parameters. 
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Figure 1 1 : Percent remaining fatigue life histogram after developing a lead crack. Case 3: 
variable unsymmetric crack sizes; mean crack growth parameters for each panel. 



Figure 12: Percent remaining fatigue life histogram after developing a lead crack, Case 4: 
variable symmetric crack sizes; variable crack growth parameters for each panel. 
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Figure 13: Percent remaining fatigue life histogram after developing a lead crack, Case 5: 

variable unsymmetric crack sizes; variable crack growth parameters for each panel. 
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Figure 14: Percent remaining fatigue life histogram. Case 6: variable symmetric crack sizes; 

variable crack growth parameters at each hole for each panel (heterogeneous panel). 
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Figure 16 : Analytically derived initial crack size distribution compared with observed cracks 
reported by Lincoln. 
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